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Interest in the elastic properties of regular lattices constructed from domain walls has recently been 
motivated by cosmological applications as solid dark energy. This work investigates the particularly 
simple examples of triangular, hexagonal and square lattices in two dimensions and a variety of 
more complicated lattices in three dimensions which have cubic symmetry. The relevant rigidity 
coefficients are computed taking into account non-affine perturbations where necessary, and these 
are used to evaluate the propagation velocity for any macroscopic scale perturbation mode. Using 
this information we assess the stability of the various configurations. It is found that triangular 
lattices are isotropic and stable, whereas hexagonal lattices are unstable. It is argued that the 
simple orthonormal cases of a square in two dimensions and the cube in three are stable, except 
to perturbations of infinite extent. We also find that the more complicated case of a rhombic 
dodecahedral lattice is stable, except to the existence of transverse modes in certain directions, 
whereas a lattice formed from truncated octahedra is unstable. 

PACS numbers: 


I. INTRODUCTION 

The possibility of an effective negative macroscopic pressure due to a medium constituted from a lattice of cosmic 
strings, and a fortiori domain walls, has motivated interest in such structures as the mechanism [IQ which might 
account for the observed acceleration of the universe U- A crucial issue that arises is whether such structures can 
be sufficiently rigid to offset the negative pressure, and behave as a stable elastic solid. A recent investigation 0 has 
shown that a randomly orientated, and thus macroscopically isotropic, distribution of cosmic strings or membranes 
will indeed be sufhciently rigid if only affine deformations are considered. It was explained that this should apply in 
the case where the junctions are of even-type (X-type or *-type), that is, when an even number of elements (strings 
or walls) meet at each junction. 

The purpose of the present work is to initiate a study of the macroscopic elastic properties, and hence the local 
stability properties, of more realistic models in which the solid is formed from a regular lattice structure, which is 
naturally anisotropic. Since wall lattices, with P/p = —2/3 in 3 dimensions, are favoured by observations we will 
concentrate on this possibility, but our methods could equally well be applied to the string case with P/p = —1/3. 
We will consider the cases of wall lattices with triangular, hexagonal and square symmetry in two dimensions, and 
three dimensional lattices with cubic symmetry which can be created from the Wigner-Seitz cells of simple cubic, 
face-centred cubic and body-centred cubic lattices. The unit cells of these primitive lattices are the simple cube, the 
rhombic dodecahedron and the truncated octahedron respectively. 

The important issue which we need to address is that some of the regular lattices we would like to consider contain 
junctions that are of odd-type (an odd number of elements meeting at a junction). Y-type junctions of this kind 
appear in, for example, models with hexagonal symmetry where three walls meet at 120 degrees, and are ubiquitous 
in three dimensional models. As was pointed out in ref. such models require consideration non-affine deformations 
to maintain the equilibrium condition at the junction, violating the assumptions of the earlier calculation. Typically 
such non-affine deformations will lead to a reduction in the rigidity when compared to the affine case, making the 
structures less stable. In addition to this, it will also be necessary to introduce an extra shear modulus to take into 
account the anisotropic nature of lattices with square and cubic symmetries. 

First, we will establish the general conditions for stability in the isotropic case and cubic cases in general dimensions. 
This formalism can then be applied to the cases of triangular and hexagonal symmetry in two dimensions. It is shown 
that such systems are effectively isotropic, but that the triangular case is stable, whereas the hexagonal case is not. 
We then apply this to simple orthonormal cases of a square in two dimensions and simple cubic system in three 
dimensions. We hnd that such systems are anisotropic and have longitudinal zero modes corresponding to a finite 
number of specihc directions. We argue that such structures are only unstable to perturbations of infinite extent. 
Finally, we consider wall lattices whose unit cells are the truncated octahedron and the rhombic dodecahedron. The 
former is unstable, but the latter is only unstable to transverse zero modes in certain directions. 
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II. FORMALISM 


A. Quasi Hookean ansatz 


The formulation of a perfectly elastic solid model is based on a material base space which is the quotient of the 
spacetime by the congruence of the material worldlines. The D-dimensional, time dependent metric induced on this 
base space is denoted 706 . For the purooses of our discussion it will be sufficient to use models which are characterised 
by a quasi Hookean equation of state [3 • This implies that the energy density can be expressed in terms of the rigidity 
tensor, 

p = ( 1 ) 


The overhead check symbol is used to indicate quantities which depend only on the value of the conserved number 
density, n. The quantity p is the minimum value taken by the mass-energy density, p, for the equivalent value n. It 
is important to distinguish between ■jab and its value jab in the reference state for which this minimum is attained. 
A quadratic form, such as o, can be expected to be valid as a realistic approximation for sufficiently small values of 
the constant volume shear tensor Sab which is defined by 

^ab — 2 ("Tofe bab) • (^) 

The functional dependence of p on n determines corresponding reference state value of the pressure tensor, 
and also of the bulk modulus, (3. These are given by 


pab ^ 


- —1 ab 


P = n- - p, 

an 


$ = n 


dP 

dn 


( 3 ) 


in terms of the inverse metric 7 induced on the base space. The rigidity tensor on the material base space has 
the following symmetry conditions 


^abcd _ ^{ab){cd) _ ^cdab 


( 4 ) 


The fact that jab and jab must have the same determinant (due to the conservation of volume) implies that on 
the D-dimensional material base the symmetric tensor Sab will have only (D^ -|- D — 2)/2 (instead of D(D -|- l)/2) 
independent components, and more specihcally that to first order it will be trace free with respect to either the 
actual metric jab or the reference metric jab- Therefore, it is necessary to impose a restriction to completely fix the 
specification of This can be most conveniently achieved by requiring that it be trace free with respect to 

the reference metric 


= 0 . 


( 5 ) 


B. Characteristic propagation equation 


We shall be considering systems characterised by the absence of compressional distortion which means that the 
only effect of a change in n on the reference metric jab is to multiply it by a simple conformal factor. Under these 
circumstances the complete elasticity tensor in the reference state will differ from the rigidity tensor by 

an isotropic contribution determined just by the dependence of the reference density p on n. This is given by 

pabcd ^ ^abcd 2pj-^ c(a^-l b)d cd ^ ('0) 


This leads to an expression of the form 

jpbcd _ ^abcd _|_ ab^-1 cd _ pp-1 a[6^-l c]d ^ 

for the relativistic Hadamard elasticity tensor 0 that is needed for the formulation of the characteristic equation 
governing the propagation of small perturbations. It follows that the characteristic equation for the speed u of a 
mode with polarisation direction La propagating in the direction specified by a unit space vector will be given by 


u 2 (p + P)7-1“'’ 



( 8 ) 
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where 


gab ^ (g) 

The eigenmodes of this equation will yield the normal modes of the system and will allow us to establish stability 
conditions. 


C. Isotropic case 


A system of the kind we have been considering will evidently reduce to one of simple perfect fluid type if and only 
if the rigidity tensor is zero, a condition that can be seen to be equivalent to that of vanishing of the scalar 

rigidity modulus fi as defined by the formula 




1 


2{D + 2){D-l) 




( 10 ) 


The normalisation in m is determined by the requirement of consistency with traditional terminology 0 in the 
ordinary isotropic case, for which usual convention is that the rigidity term in m should simply be given by 

1 ^abcd - — lc(a - lb)d ^ ab 

-E SabScd = ^l ^ 7 SabScd = fJ-SabS , ( 11 ) 

so that the corresponding expression for the isotropic rigidity tensor will be 

^abcd ^ 2fi ( 12 ) 

This means that the non vanishing Cartesian components of the rigidity tensor in 0 will just be of three 

kinds. There will be the mixed kind specified by a pair of axes, Ai ^ ^ ^, the crossed kind also specified by a pair of axes, 
and finally the pure kind specified by a single axis, According to m the values of these components 

will be given by 




E 


112 2 _ 


2jl 

')D ’ 


17““ = 


2{D-l)fL 

D 


(13) 


It can be seen from 113 that, in a medium of this isotropic kind, the polarization eigenvector La in 0 will have one 
of two distinct types. It may be orthogonal to the propagation direction i/“, iaV°' = 0, in which case the corresponding 
transverse velocity eigenvalue ut will be given simply by 


2 _ A 

'Urp - 


P + P 


(14) 


Alternatively the polarisation direction ta may be parallel to the propagation direction in which case the corre¬ 
sponding longitudinal velocity eigenvalue ul will be given by 


^^^Dp + 2{D-l)fi 


D{p + P) 


(15) 


There are D — 1 orthogonal directions in the eigenspace with velocity vt and one with ul. In the case where the 
equation of state is given by P = wp the two sound speeds are given by 


'V'-p — 




(1 -I- w)p ’ 


Vt =w + 


2{D-l)fi 
D{1 + w)p ' 


(16) 


We shall be concerned with cases where the medium will consist of a relativistic system of Dirac-Goto-Nambu type 
branes. In such a system the reference state can be described by a polytropic equation of state with index 7 , so that 


P = wp , (3 = "/wp , w = 7 — 1, 


(17) 
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where 7 is the ratio of the codimension to the space dimension, so that for the generic p-brane case one obtains 

(D-p) 


7 = 


D 


P 

w = --. 


(18) 


This means that both for strings, as characterised by p = 1, and walls, as characterised by p = U — 1 we shall always 
have 


{D-l)p 


£)2 


(19) 


for arbitrary space dimension D. Thus, in the ordinary 3 dimensional case the bulk modulus will be given by 
(3 = —2pl9, not only for strings, with 7 = 2/3 and w = —1/3, but also for walls, with 7 = 1/3 and w = —2/3. In the 
2 dimensional case, for which there is no need to distinguish between strings and walls, one will simply have /3 = —p/d 
with 7 = 1/2 and w = — 1 / 2 . 

In all these cases, it can be seen that for a system of the isotropic kind characterised by the corresponding 
propagation speeds of transverse and longitudinally polarised perturbation modes will be given by 



, 2{D-1) fjl 1 \ 

7 ^ U ‘^Dj ■ 


( 20 ) 


D. Stability and causality in the isotropic case 

In order for the system to be stable and respect causality, we require that the two sound speeds satisfy 0 < < I 

and 0 < < I, since otherwise there will be modes propagating at speeds faster than the speed of light, or there will 

be modes whose amplitude grows exponentially on a timescale proportional to its wavelength. In the general case of 
P = wp this implies that 


0 < - < I 


w, 


Dw{l + w) ^ fi ^ D{1 — w'^) 
2{D-l) -J)- 2{D-1) ’ 


( 21 ) 


where we have assumed that 1 + w > 0 . 

For systems described by (I17ll8l19l2()li . the stability criterion can be seen to be that the rigidity modulus is not 
less than a critical minimum value given 



( 22 ) 


which implies that the critical minimum value is /t = p /6 for D = 3 and /t = p/4 for D = 2. For an isotropically 
randomised distribution of strings or membranes considered in our previous work 0 this will be satisfied since it was 
shown that the relevant rigidity modulus will be given for D = 3 hy fi = 4p/15 , and it can easily be verified that 
an an analogous calculation for the simpler case D = 2 gives fi = 3p/8. Moreover, the condition for causal wave 
propagation is that 



A ^ 1 yD 
J) ^ 2{D - 1) ■ 


(23) 


In cases of present interest with D = 2 and D = 3, it can be seen that the first of these two formula yields the 
strongest constraint, but this need not be the case for larger values of H. It is clear that these are also satisfied by 
the model discussed in ref. 0 

We should note that we have allowed for the possibility of zero modes in our analysis above, that is, we have classed 
the cases of wt = 0 and wl = 0 as stable. This is something which is somewhat delicate and more consideration 
of that particular context is necessary. A familiar example of such a mode is provided by the ordinary barotropic 
perfect fluid, as characterized by a strictly positive bulk modulus, /3, and vanishing rigidity modulus, p. In this case 
> 0, but vt = 0. This indicates that transverse shear modes do not propagate and are not subject to exponential 
growth. In a non-expanding background, they are subject to unbounded linear growth and therefore can be described 
as marginally unstable. However, in an expanding background they are expected to be subjected to either power law 
growth or decay dependent on the details of the situation. 

Furthermore, it is worthwhile to reflect on how the existence, in principle, of such marginally unstable modes can 
be compatible with our everyday experience of the stability of ordinary static fluid states. The key to the paradox is 
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the realisation that the modes in question are plane waves of infinite extent. For such an instability to grow, there 
must exist growing perturbations that are initially confined within a locally bounded support region. The existence 
of marginally unstable modes of plane wave type and thus with unbounded support is not always sufficient to provide 
unstable wave-packet combinations with bounded support. 

An example in which there are no exponentially unstable modes, but in which (unlike the perfect fluid case) there 
actually are sufficient marginally unstable plane wave type modes to form growing perturbations with local support 
is provided by the hexagonal lattice model discussed in the next section. This Y-type example is to be contrasted 
with the X-type cases discussed in subsequent sections, namely the square and cubic lattice models, for which some 
zero modes do exist but, as in the fluid case, are not enough to give rise to an effective locally supported instability. 


E. Simple anisotropic cases 


The possible anisotropic forms for the rigidity tensor can be classified for Z? = 2 and Z? = 3 in terms of the 
Bravais lattices @ extensively studied in crystallography. We will consider here the simplest of these cases, that of 
the orthonormal type, meaning lattices of simple square type for D = 2 and lattices of simple cubic type for ZZ = 3, 
not to mention hypercubic types for larger values of D that might be worth considering in view of the current interest 
in higher dimensional cosmological models. It is well documented in the literature that, due to their symmetry under 
permutation of the axes, models of this type have non vanishing Cartesian components of the rigidity tensor 
in o of three kinds, with just two independent rigidity parameters. We will suggestively designate these parameters 
as p-T and /II- There will be the mixed kind specified by a pair of axes, for example the crossed kind also 

specified by a pair of axes, for example and finally the pure kind specified by a single axis, for example 

The requirement 0 that the rigidity tensor be trace free implies that the value of a component of the pure type must 
be given in terms of that of a component of the crossed type by a relation of the form 


y;iiii = (1 _ ZI)r“ = " . 


(24) 


In order to ensure that they both go over to the ordinary scalar rigidity modulus p in the isotropic limit as given by 
dnj the choice of the independent rigidity coefficients pT and pi^ can thus be conveniently fixed by the specifications 




= ^J-T : 




2Al 
' D 




2(ZZ-1)Al 

D 


(25) 


which entail that the rigidity scalar p defined by will be given, as a weighted root square mean, by 


~2 _ PpT + 2 / 2 ^ 
^ ~ D + 2 


(26) 


Moreover, we define Ap = p^ — pi^ which quantifies the anisotropy of the lattice. 

Writing = cos9a, where 9a is the angle between the propagation direction and the x® axis, it can be seen from 
the formula 0 that the diagonal components of the matrix in the characteristic equation © can be written as 
(no summation) 


= 


D 


pL+ P] COS^ 9i + pT sin^ 9i, 


for 1 < * < ZI and the off-diagonal components are given by 


_ j --jy + At + /3 ) cos 9i cos 9j , 


(27) 


(28) 


for j. 


F. Stability and causality in the general orthonormal case 

In an anisotropic system, such as the cubic type considered here, the propagation speed of perturbations modes 
are typically spatially dependent making the discussion of stability more complicated than in the isotropic case. The 
characteristic propagation equation can be solved by solving an eigenvalue problem for \ = {p + P)v^. In general 
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this is a complicated polynomial equation of degree D, with the coefficients depending on cos0i for i = The 

characteristic equation for general D is complicated; here we present it for the cases of specific interest D = 2, 

3A^ — (4/iL + 3/3 + 3/iT)A + + 4/fL) + 4i32(/iL ~ /fT)(AL + 3/3) = 0, (29) 


and 13 = 3, 


3A^ — (4/fL + 3/3 + 6 /iT)A^ 

+ {/tT(8/iL + 6/3 + 3/iT) + 4(/iL - AT)(/iL + 3 / 3 ) 53 } A (30) 

~I^T (4Al + 3/3) — 4/iT(/fL — /3t)(Al + 3/3)53 — 4(/iL — /iT)^(/iT + 3/3)C = 0, 

where B 2 = cos^ cos^ O 2 , = cos^ 9i cos^ 82 + cos^ 9i cos^ 6*3 + cos^ cos^ 6*3 and C = cos^ 9i cos^ 6*2 cos^ 6*3 

For specific value of the rigidity coefficients /11 and /fx the D solutions of the characteristic equation can be 
computed and unlike the isotropic case these will, in general, be distinct and be dependent on direction. 

We have been able to solve (P|l in general, 

2 (4/iL + 3/3 + 3/iT) ± \J (4/iL + 3/3 — — 4852 (/11 — /iT)(AL + 3/3) 

which yields stability conditions /jl > /5/4 and /ix > /i/4 for /3 = —/i/4 as suggested by (IT!^ . However, in the case 
5 = 3 and, indeed, for larger values of D we have been unable to solve the characteristic equation analytically. Of 
interest when considering the stability are the regions of the (/iL,/ix) space where A = 0 is a solution. This will 
require the coefficient of A® in the characteristic equation to be zero. For general 5, this coefficient is minimized for 
two specific directions: (1) D axial directions corresponding to 0^ = 0 for each 1 < / < 5 and 9j = t: 12 for i ^ j, and 
(2) a single non-axial direction cos0i = 1/5 for all i = 1..5. 

The characteristic equation for general 5 can then be factorized. In the axial directions the characteristic equation 
is given by 


/ix — A 


£ 1-1 


2(5 - l)/iL + 5/3 - 5A =0 


and in the non-axial direction it is given by 


(32) 


2/iL + {D — 2)/ix — 2A 


D-l 


2(5- l)/ix-b5/3-2A =0 


(33) 


If we assume that $ is given by then one can ensure that all propagation modes are stable in these directions, 
if the rigidity moduli are not be lower than critical values given by 


Al = At = ^, (34) 

which is reminiscent of the isotropic case and is exactly the same as for 5 = 2. 

Therefore, we have shown that if the rigidity moduli are given by El then there are non-propagating modes. 
We have not shown rigorously that these are the unique points at which there is a zero eigenvalue solution to the 
characteristic equation, nor have we shown that that there is an negative eigenvalue, and hence an unstable mode 
below this point. However, we have checked this numerically by computing the lowest eigenvalue for a discrete set of 
points in the (/iL,//x) plane, confirming that the conditions for stability are that 


^ > J_ ^ 

p - 2D’ p - 2D' 

Similarly, one can show that causality requires that p.'i,/p < 7 and pt/ p < 7 for 5 = 2 and 5 = 3. 


(35) 


III. SHEAR MODULI FOR TRIANGULAR AND HEXAGONAL LATTICES IN 2D 

A particularly instructive - simple, but non-trivial - example is provided by the case of a regular lattice of Dirac- 
Nambu-Goto type domain wall boundaries separating hexagonal domains in 2 dimensions. It is illuminating to 
compare it with the more artificial and technically more trivial case of its dual lattice, the simple triangular lattice. 
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FIG. 1: The case of a triangular cosmic string lattice with even type “asterisk” form junctions in 2 dimensions, a system that 
is quite stable, as it admits no localised perturbations with negative or zero energy. 


In both the triangular case, as illustrated in Figure 1, and the hexagonal case, as illustrated in Figure 2, the 
length and direction of each boundary section is represented by one, or other, of three distinct 2 dimensional vectors, 
with components 6“^, 6“ which branch off from each other at 120° angles. Without loss of generality the space 

coordinates can be chosen in such a way that in the unperturbed state the first of these is a unit vector negatively 
oriented along the axis, so that the complete set will be given by 

!>,; «(- 1 . 0 }. ( 36 ) 


It is evident that such boundary segments might be fitted together to form a regular periodic triangular lattice, with 
any two of the three boundary segment vectors as a fundamental dyadic pair of independent symmetry generators, 
and with junctions of the, six fold, even “asterisk” type (*-type). Such a junction between an even number, in this 
case 6, of terminating segments can also be described as the crossing of half that number, in this case 3, uninterrupted 
boundaries. It follows that its equilibrium will be preserved by a simple affine linear space deformation. 

Although this network appears to be anisotropic, it is in fact not the case. To see this, let us evaluate which are 
the non-vanishing components of the rigidity tensor. Since the system has an invariance with respect to rotations, 
the easiest way is to look at the problem in the complex coordinates as defined by 

^ = x + iy^ T] = x-iy, (37) 

where (x,y) are the Cartesian coordinates. In these new coordinates, a rotation of angle a comes back to the simple 
transformation 




rj ^ e *“77 


(38) 


Since the system is invariant under rotations of angle a = ^ the only non vanishing components of the rigidity 
tensor will be the ones with as many ^ as 77 , that is and . However, the fact that is traceless brings 

back the non vanishing components to one, which is easily seen to be . This implies that the system is isotropic. 
It is to be remarked that the previous argument applies as well to the Y-type junction hexagonal 2D lattice (which is 
invariant under rotations of a = ^ ). By evaluating the effect of the affine deformation on the triangular lattice in a 
similar way to our preceding work , it can be shown that the rigidity modulus is given by the same formula as for 
an isotropic random distribution discussed earlier, namely 



(39) 


The fact that this rigidity value is substantially larger than the relevant two dimensional critical value, jl = p/4 ensures 
that such an even type triangular lattice will be absolutely stable. It can be seen from that the corresponding 
propagation speeds for transverse and longitudinally polarised perturbation modes are given by 


V 


2 _ 

T — 


3 

4 ’ 



(40) 


FIG. 2: The dotted arrows 
a regular hexagonal lattice 


indicate how the lattice periodicity generators ^ “ are placed in relation to the material structure of 
of domain walls with even type Y-form junctions in 2 dimensions 


While this stability feature may be considered to be an advantage from an engineering point of view, on the other 
hand from a more purely physical point of view, such a triangular lattice has the drawback of being relatively difficult 
to obtain in a natural context. This contrasts with the case of a hexagonal lattice such as can be formed naturally 
by a system with three symmetrically related vacuum states. 

A regular hexagonal lattice can be constructed by fitting together boundary segments of the same three kinds as 
in the triangular case. However, instead of being directly identified with any two individual members of the set three 
boundary segment vectors 6“^, 6“ as in the triangular case, the fundamental dyad of periodicity generating lattice 

vectors, and £“ say, can be taken to be any two differences between pairs of such boundary segment vectors. For 
example, one can take 



( 41 ) 


The fact that such a lattice is of odd type, with triply intersecting Y-type junctions, means that when it is subject 
to a static macroscopically affine deformation, as specified by a linear transformation of the form 

e2-^il = {S^b+e%)i2, ( 42 ) 

acting on the lattice periodicity generators, it is no longer possible to suppose that this affine transformation acts 
locally on the mesoscopic structure characterised by the boundary segment vectors 6“^, Thus, in order to find 

the new equilibrium state, which will allow us to compute the rigidity modulus, we need to find the transformations 
that, when applied to the boundary vectors ensure that all the periodicity vectors given by (EH, are subject to the 
same affine transformation (1^ . Although this is a difficult problem in general, there is a simple solution in the 
particular case of the hexagonal lattice. This involves changing the lengths, but not the directions, of these boundary 
segment vectors, according to 


% ^ 0 ) - (1 + ^( 0 ))^( 0 ) > b2 - {1 + ejb2 ■ ( 43 ) 

These transformations, which obviously preserve the equilibrium condition, will give rise to a fractional density change 
that will be given by 


^ _ £(o) +£++£_ 

P 3 


( 44 ) 


The effect of the non-affine transformation 63 on the lattice vectors specified by (EH can be expressed as an affine 
transformation of the form (1^ . where the deformation matrix components are given by 


4:6 


= 


(0) 


+ s_ 


6 ^ = 


+ S_ 


2 


2^3 


6 


£■ 


2 


= 


( 45 ) 
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In order to evaluate the components of the rigidity tensor we only need to consider transformations that conserve 
volume, a condition which can be expressed by the requirement that the determinant of the affine transformation 
matrix should have unit value. This implies the restriction 


£( 0 ) +£++£- 


1 

- 2 '" 


(o)£+ +£(o)£- +£+£-) 


(46) 


This condition enables us to eliminate and work just with the two remaining independent parameters and e_, 
thereby expressing the deformation matrix, to first order in the deformation amplitude, in the form 

e% = s\ + 0{e^}. (47) 


Here, is the trace free matrix with components given by 

£. T £ 

gi =-s2 =--^ 


e,—e 
= sb = — - 


(48) 


2 ’ “ 1 “ 2 2^3 

In terms of these shear matrix components it can be seen that the fractional density variation llH) can be written as 


which is manifestly of the isotropic form m with a rigidity coefficient that can be deduce to be 


(49) 



(50) 


In the discussion above, we chose transformations that would preserve equilibrium and would be consistent with 
the periodicity vectors being subject to the same affine deformation, that is, working from the mesoscopic to the 
macroscopic scale — a bottom-to-top approach. We should point out that we could also have done the calculation 
the other way around in a top-bottom approach. In this case, the problem is posed in the following way: consider a 
Y-junction and subject all of its points (excluding the vertex) to a given deformation and then look for the position 
of the vertex which, keeping the new endpoints of the boundary vectors fixed, minimizes the energy. If we consider 
the origin and the points b'^ as defining the Y-type junction then a transformation of the 6“s under x —*■ (1 + e)x 
and y y/{l + e) and the origin moving to the point (x, y) results in an energy functional 


E{x,y) 



V3 

2 ( 1 + £) 




V3 

2 ( 1 + £) 



(51) 


One can then compute the position (x,?/) which minimizes E{x,y) either numerically or, since we only need the 
positions to order i5^, as a power series. We find that x = e -I- e^/2 -|- 0{£^} and y = 0{£^} are the minima and so 
Ap/p =1/2 and we conclude the that fi/p =1/4. One can also consider a simple shear transformation x ^ x -|- 2e?/, 
y ^ y which yields 



-k x/oTxpT^, 


(52) 


which is minimized by x = 2£^ + 0{£^} and y = £ + 0{£^} and yields the same result. While the bottom-top approach 
is elegant, it is much more complicated when dealing with 3 dimensional lattices and is often intractable. As we shall 
see in section m in most cases the top-to-bottom approach is more practical, but even then it is often necessary to 
resort to numerical means to solve for the minimum point. 

It can be seen from (I19II that /3 = —p/4 for a, D = 2 string/wall system and therefore using (1151) the value isni) 
obtained for the rigidity modulus of this odd type hexagonal lattice has exactly the critical value for stability of 
longitudinal modes. According to da and m the squared velocities of the transverse and longitudinal modes will 
be given respectively by 


V 


2 

T 


1 

2 ’ 


vl=0. 


(53) 
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FIG. 3: Locally confined energy conserving, and therefore linearly destabilising, deformation in “odd” type hexagonal lattice 
of domain walls with Y-type junctions in 2 dimensions. 



FIG. 4: Innocuous " albeit energy conserving “ deformation, of axially aligned plane wave plane wave type, in square lattice 
with even type X form junctions in 2 dimensions. Unlike the example in Figure El such a perturbation is not locally confined, 
since it can be seen to be incompatible with fixed boundary conditions at the upper and lower borders of the figure, so it does 
not engender any local instability. 


This contrasts with the corresponding case iBUI) of the even type triangular lattice for which, the value of the 
rigidity modulus was found to be 50 percent larger, and so safely within the range required for absolute stability. 

It is instructive to compare this hexagonal lattice case, for which the marginal zero velocity modes are of (irro- 
tational) longitudinally polarised type, with the case of a simple perfect fluid for which, as remarked above, the 
marginal zero velocity modes are of transversely polarised (density conserving) type. The latter are not sufficient for 
local destabilisation of the fluid. However, it can be seen that the hexagonal lattice model really will be linearly, but 
not exponentially destabilised by its zero velocity longitudinal modes since they exist for all wavenumber directions, 
they can form a sufficient basis for superposition to constitute irrotational wave packets with locally confined support. 
The mesoscopic structure of an example of such a locally confined, zero energy, and therefore linearly destabilising, 
perturbation in a regular hexagonal lattice is shown in FigureEl This shows the surprising result that although duality 
preserves the isotropic behaviour, it is not so for the stability. 


IV. ORTHONORMAL LATTICES 


We now turn our attention to orthonormal lattices which vt\ D = 2 correspond to lattices with square symmetry and 
in U = 3 to those with cubic symmetry. The analysis of such cases is facilitated by the fact that their junctions are of 
even type - more specifically their walls meet at X type crossovers - which means the effect of a static macroscopically 
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affine deformation can be represented by an affine transformation even at a local mesoscopic level. In particular this 
is because the X-type junctions are stable for any opening angle, as opposed to the very specific angle (120°) in the 
hexagonal case. 

To actually evaluate the relevant Cartesian rigidity components for the orthonormal lattices, we consider the effect 
of affine transformations of the general form 

a;® 1-^ X* = x* + s '', (54) 


as specified by a deformation matrix whose symmetrized trace free part. 


Sj + e-® 






D 


(55) 


can be interpreted as the Cartesian representation of the shear tensor given by ©• 

To evaluate the mixed component, consider the volume conserving simple shear deformation for which 

the only non vanishing component of the deformation matrix is e® = 2 e say, so that the only non vanishing shear 
components of the corresponding shear matrix will be = e. This transformation has no effect on any of 

the relevant (originally) orthonormal basis vectors except the first, whose length will change by an amount 5^ = 2e^, 
thereby bringing about a fractional change 5p/p = 2e^ jD in the mean density. Since the effect of the corresponding 
shear deformation llK^ . whose only non-vanishing matrix components are =5^1 = £ will be given according to o 
by &p = 2X1 = ®^ , the required shear component value can be deduced to be 


^1212 


P_ 

D ■ 


(56) 


To evaluate the other independent component, consider the irrotational volume conserving shear deformation for 
which the only non vanishing components of the deformation matrix are given to quadratic order by = e, e ^2 — 
—e + e^. In this case there are two (originally) orthonormal basis vectors whose length is affected, namely the first, 
for which Si = e, and the second, for which Si = —e + The combined effect of these length adjustments evidently 
cancels out to first order, producing a fractional density change that is given at quadratic order by Sp/p = e'^fD. 
In this case the only non-vanishing shear matrix components will be given to first order by s\ = —s\ = e, so the 
corresponding density change will be given according to (j^J by Sp/p = (27“ “ — 27®“^)e^. It follows from (1241) that 
the required values can be read out as 


27““ = 


jP - 1)P 
£)2 


27 


112 2 _ 


P 

'D 2 • 


(57) 


The outcome, therefore, is that for a lattice of the orthonormal kind, the rigidity coefficients introduced in m will 
be given by 


- P - P 

Using our earlier analysis, it is clear that the longitudinal rigidity is only sufficient for those modes to be a marginally 
stable. We can also deduce, the scalar shear modulus 


P_ 2D + 1 
~ V 2 D^{D + 2) ’ 


(59) 


and the anisotropy parameter is 


A/2 1 


(60) 


For D = 2, we have that (i/p^ 0.395 and Ap/p « 0.250, whereas for Z7 = 3 one can deduce that p/p « 0.279 and 
A/i/p « 0.167. 

When used in conjunction with the formula (tT^ for the bulk modulus, it can be seen that the relations (CTl lead 
to a dramatic simplification of the matrix Q°'^ in the characteristic equations as given by the of the formulae (| 77 |) 
and in which most of the terms cancel out, leaving 




D 


sin^0. 


Q'^ = 0 . 


(61) 
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This implies that the matrix Q°'^ will be of purely diagonal form, which implies that, unlike the isotropic case, the 
allowed directions of the polarisation eigenvector in are independent of the propagation direction, being fixed 
to coincide with the directions of the preferred orthonormal axes. On the other hand, still unlike what occurs in 
isotropic cases, the eigenvalues and hence the propagation speed will depend on its direction. 

In particular, by substituting ED into it can be seen, for example, that the mode with aligned with the x* 
axis will propagate with a value say that will be given by 


V 


2 

(0 


sin^0. 

'fD 


(62) 


for i = 1 ,.., £) where 7 is the relevant polytropic index, which according to CD will be given by 7 = 1/D for a lattice 
of cosmic domain walls and hy j = (D — 1)/D for the other case to which the present analysis applies, namely a 
lattice of cosmic strings. 

In the physically interesting case of domain walls, independent of the dimension we end up with the memorable 
formula 


= sin^ 6 »i, (63) 

which tells us that the velocity ranges up to the speed of light in the limit of propagation orthogonal to the polarisation 
direction, and down to zero for propagation parallel to the polarisation direction. In the more artificial case of an 
orthonormal string lattice the velocity will also range down to zero for propagation parallel to the polarisation direction, 
but (except in the 2 dimensional case where the distinction between walls and strings does not arise) the maximum 
velocity possible for propagation orthogonal to the polarisation direction will be reduced by a factor l/^/D — 1. 

Unlike the highly stable example of the triangular lattice, which according to (EUll has no zero velocity modes at all, 
the orthonormal lattices considered here do have some exceptional zero velocity unstable plane wave modes, namely 
those with an exactly vanishing value of the angle in measuring the deviation of the propagation from the relevant 
preferred axis. However, as all the neighbouring modes are strictly stable, the occurrence of these exceptional modes 
is not sufficient for the construction of unstable locally confined perturbations. The situation in these square and 
cubic examples is therefore similar to that in an ordinary fluid, which has unconfined linearly unstable shear modes, 
but is nevertheless effectively stable with respect to perturbations having confined support. These modes correspond 
to the kind of global deformation illustrated in Fig. 0| 

This contrasts with the effectively unstable case of the hexagonal lattice for which mesoscopically non-afline devia¬ 
tions reduce the energy associated with a macroscopically affine deviation to such an extent that, according to (I53II it 
possesses linearly unstable plane wave modes of longitudinally polarised type in all directions. Although such waves 
are of individually unconfined support, they are sufficiently numerous for their collective superpositions to provide 
linearly unstable perturbations that are initially confined within a finite region, thereby effectively destroying the 
stability of the system as a whole. 


V. PRIMITIVE LATTICES IN 3D WITH CUBIC SYMMETRY 

We have already commented that the anisotropic forms of the elasticity tensor can be classified in 3D by the Bravais 
lattices. Since our ultimate aim is the application to the idea of Solid Dark Energy, it seems sensible to consider the 
three primitive lattices allowed in 3D which are the Wigner-Seitz cells of simple cubic (SC), body-centred cubic (BCC) 
and face-centred cubic (FCC) lattices. The polyhedra corresponding to these lattices are the simple cube, already 
studied in the previous section, the truncated octahedron (sometimes known as the the tetrakaidecahedron |l(T|'l and 
the rhombic dodecahedron. The latter two fall into the class of Archimedean solids and are space-filling, as is the 
former which is a platonic solid. The truncated octahedron has 14 faces, 6 of which are squares and 8 are hexagons. 
The rhombic dodecahedra is a polyhedron with 12 faces each of which is a rhombus. 

Unfortunately, the analytic techniques applied in the previous sections are intractable in these cases. However, since 
similar ideas have been studied extensively in the context of soap films, foams and other soft condensed matter |ll|. a 
software package, known as the Surface Evolver 0 , has been developed to minimize the surface area of polyhedra of 
arbitrary topology. This package, which is publicly available, allows one to accurately compute the energy of minimal 
energy configurations in an unperturbed state and ones which are subject to a small deformation. 

To obtain the rigidity coefficients we select deformations in an analogous fashion to the previous sections. The mixed 
component of the rigidity tensor is evaluated by considering the volume conserving simple shear deformation 

for which the only non vanishing component of the deformation matrix is = 2e; we will denote this DEFl in the 
subsequent discussion. This will allow us to compute /xx. To evaluate the other independent component, /xl, we 
consider the irrotational volume conserving shear deformation for which the only non vanishing components of the 
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FIG. 5: A lattice of rhombic dodecahedra in an unperturbed minimal energy configuration. This configuration with fiat faces 
is the minimal energy configuration. 


deformation matrix are given to quadratic order by e\ = e, e ^2 = —£ + this will be denoted DEF2. In both cases 
we compute the fractional energy change Sp/p for a range of small values of e and fit for the coefficient of e^. 

Let us first consider a lattice of rhombic dodecahedra in an FCC arrangement as illustrated in FigEl The config¬ 
uration has flat faces and by using Surface Fvolver we find that this is the local minimum energy configuration with 
this topology. This structure contains a significant fraction of odd-type junctions as well as some even type junctions 
which will help to stabilise the system as they will not be required to adjust themselves after an affine transformation 
- in a system with purely even-type junctions the non-affine correction is zero. 

Allowing for only affine deformations with of the kind DFFl, we find that the fractional change in the mean density 
is 6 p/p = O.Se^. Subsequent energy minimisation using the Surface Fvolver gives Sp/p = O.SSSe:^. Using DFF2, we 
find that Sp/p = 0.583e^ and Sp/p = 0.445e^ for the affine and non-affine cases, respectively. Using these values we 
find that 


— « 0.222, 0.167, (64) 

P P 

for this configuration when we take into account non-affine perturbations. Using higher numerical accuracy than the 
3 decimal places quoted here, we deduce that it is likely that pi^/p = 2/9 and pt/ p =1/6 although this can only be 
checked to a particular level of accuracy. 

Next we consider a lattice of truncated octahedra in BCC arrangement. This is the configuration which was 
proposed by Lord Kelvin as the polyhedron with the lowest surface area for a fixed volume — the so called Kelvin 
Problem 0,^3- Although this is now known not to be the absolute minimum M , it is has been extensively studied 
in the context of soap films. As illustrated in Fig. 1^1 truncated octahedra with flat faces are not the local minimum 
energy configuration: the minimum has curved faces. 

Under an affine transformation DEFl there is a fractional change in the mean density of Sp/p = 0.494e^ for 
truncated octahedra with flat faces. A subsequent energy minimisation using the Surface Evolver allowing for non- 
affine transformations, reduces the fractional change in the mean density which is given by Sp/p = 0.367e^. For 
DEF2, the fractional change in the mean density is Sp/p = 0.592£^ for the affine deformation - subsequent vertex 
readjustment due to energy minimisation reduces this to Sp/p = 0.218e^. Therefore, we find that for the Kelvin foam 
of truncated octahedra the rigidity coefficients are given by 

0.183, 0.109. (65) 

P P 

These values are compatible with previous calculations done in the context of soap films M- 

We have tabulated the compute shear moduli for the SC, FCC and BCC cells in table 13 along with the computed 
scalar shear modulus p//> and the anisotropy A/r/p. We see that the value of p/p is lower than the value computed 
in ref. 0 due to the effects of non-afRne deformation and also the strong anisotropy present in all cases. Using the 
results derived in section III FI one can see that the SC and FCC cell structures are conditionally stable with zero 
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FIG. 6: On the left is an illustration of the BCC arrangement of truncated octahedra in unperturbed state and on the right is 
the unperturbed minimal energy conhguration created using the Surface Evolver. Notice the slight curvature of the faces on 
the right. 



Ml/p 

Mt/p 

m/p 

Ap/p 

SC 

FCC 

BCC 

1/6 

2/9 

0.109 

1/3 

1/6 

0.183 

v/7/90 

v/59/1620 

0.150 

1/6 

-1/18 

0.074 


TABLE I: The shear moduli computed using for the simple cubic (SC), face-centred cubic lattices and body-centred cubic 
(BCC) lattices. 


modes in each case. These are longitudinal in the SC cases and transverse in the the FCC case. Moreover we find 
that the BCC lattice is unstable since /tl = 0.109 < 1/6. 

It is not possible to compute the sound speeds analytically for the FCC and BCC cases as we did for the SC case 
in the previous section. In Fig. [ 7 | we have plotted the numerically computed square sound speeds (v^) as a function 
of (0, (/>) for all three lattices in an equal area Hammer-Aitoff projection where k = fc(sin6icos(/), sin6isin0,cos6i). In 
all three cases the characteristic four-fold symmetry is evident. However, there are some interesting differences: for 
the BCC lattice there is a substantial region of the space (0, </) for which < 0, confirming the instability already 
identihed; the SC lattice has = 0 at the centre of each of the faces of a cube, whereas the FCC lattice has = 0 
in directions exemplified by 


VI. DISCUSSION AND CONCLUSIONS 

In this paper we have examined the stability of a variety of different lattices both in 2 and 3 dimensions (some of 
our results even apply in arbitrary dimension). There are two features which we have investigated which go beyond 
our earlier analysis [^. Firstly, we have conhrmed the physical importance of the distinction between odd type and 
even type junctions between the relevant wall (or string) segments. It appears that the more odd type junctions a 
given lattice cell has the less stable it is likely to be. The other feature which we have highlighted are the effects of 
anisotropy. Even though the orthonormal lattices contain only even type junctions, the strong anisotropy present in 
these models is important. 

Since our main aim was to find a lattice which is stable as a candidate for solid dark energy it is worth discussing 
progress toward this. The 3 dimensional configurations which we have investigated are formed from cells based around 
the simple cube, rhombic dodecahedra and the truncated icosahedron. These are the primitive Bravais lattices with 
cubic symmetry. 

Our calculations suggest that we should immediately eliminate the cell based on the truncated octahedron which 
is exponentially unstable for a wide range of (0, (/>). This is something of a surprise since it is known at least play a 
role in the structure of soap films. However, this highlights an important difference between that field and the one at 
issue here. Namely the bubbles in soap hlms have an extra contribution to the bulk modulus due to the air pressure 
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FIG. 7; The computed square sound speeds for the three orthogonal modes plotted as a function of the polar angles 
of Fourier space in the equal-area Hammer Aitoff project. The top three are for the simple cube, the middle three are for 
the lattice of rhombic dodecahedra and the bottom three are for the lattice of rhombic dodecahedra. Black on the colour bar 
denoted points where the square sound speeds are negative corresponding to exponentially unstable regions 


in the bubble, whereas we require that the domain wall lattices exists in vacuum. 

We have shown that the rhombic dodecahedral lattice exhibits transverse zero modes. The question of whether 
these zero modes are of infinite extent or have local support is more difficult as the junctions are of mixed type. As 
there are eight principle directions for which there is a zero mode we conjecture that it is possible to construct a 
closed path around these directions (as in the hexagonal case). In the hexagonal case the length scale of this finite 
mode is comparable to the elementary cell size. However, due to the existence of the even-type junctions in the 
rhombic dodecahedron the length scale of the mode may potentially be large in comparison to the elementary cell. 
We therefore conclude that this mode will be of larger extent than in the hexagonal case, although this should be 
tested numerically. Indeed, as described below, this will be the subject of some future work. 

Remarkably, the best candidate that we have at present is the simple cubic lattice which has only has unstable 
modes which are of infinite extent. These zero modes correspond to there being no specific distance fixed between 
infinite straight walls with the model based on walls with Nambu-Goto equation of state. Clearly this issue, along 
with the precise spatial distribution of the vacua, would need to be addressed by any field theoretical model. We are 
presently investigating this issue in an 0{N) model with cubic anisotropy and the results will be presented in a future 
publication M 

One should not take our comments above to reflect pessimism about the Solid Dark Energy scenario; this analysis is 
only the beginning of a search for stable configuration. There are number of possibilities for further investigation: in 
particular, the other primitive Bravais lattices with lower symmetry than cubic and compound structures with cubic 
symmetry such as the tetrahedral close packing structures. We plan to investigate these possibilities in future work. 
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